Mon. Not. R. Astron. Soc. OOP. HTOl (20011 Printed 6 October 2010 (MN WT^ style file v2.2) 



Positron annihilation spectrum from the Galactic Centre 
region observed by SPI/INTEGRAL, revisited: 
annihilation in a cooling ISM? 

. E. Churazov/'^ S. Sazonov,^ S. Tsygankov/'^ R. Sunyaev,^'^ D. Varshalovich^ 

' ^ Max- Planck- Institut fiir Astrophysik, Karl-Schwarzschild-Strasse 1, 85741 Garching, Germany 

^ Space Research Institute (IKI), Profsoyuznaya 84/32, Moscow 117997, Russia 
• ^ loffe Physical Techincal Institute, Polytekhnicheskaya 26, St Petersburg 194021, Russia 



6 October 2010 



ABSTRACT 

Wc analyse SPI/INTEGRAL data on the 511 kcV line from the Galactic Centre, 
accumulated over '^6 years of observations. We decompose the X-ray and soft gamma- 
ray emission of the central part of the Milky Way into a relatively compact "Bulge" 
and a more extended "Disk" components and report their spectral properties. The 
Bulge component shows a prominent 511 keV line and essentially no flux at 1.8 MeV, 
while the Disk component on the contrary contains a prominent 1.8 MeV line and a 
very weak annihilation line. 

We show that the spectral shape of the annihilation radiation (the narrow 511 keV 
line and the associated othro-positronium continuum) is surprisingly well described by 
a model of annihilation of hot positrons in a radiatively cooling interstellar medium 
(ISM). The model assumes that positrons are initially injected into a hot (^ 10^ K), 
volume filling ISM, which is allowed to freely cool via radiative losses. The annihilation 
time in such a medium is longer than the cooling time for temperatures higher than a 
few 10'' K. Thus, most of the positrons annihilate only after the gas has cooled down 
to ~ 10^ K, giving rise to annihilation emission characteristic of a warm, ionized ISM. 

Key words: Galaxy: centre ~ gamma rays: observations - ISM: general 



1 INTRODUCTION 

Although the annihilation line of positrons at 511 keV 
is the brightest gamma-ray line in the Galaxy, the ori- 
gin of the annihilating positrons in not yet firmly es- 
tablished. First observed with a Nal scintillator as a 
^476 keV line coming from the Galact ic Centre (GC) region 
llJohnson. Harnden. fc Havmed 1 19721 : IJohnson fc Havmed 
119731 ). it was subsequently unambiguously identified with 
a narrow (full width at half maximum, FWHM< 3.2 
keV) e^e~ annihilation line us i ng ge rmanium detectors 
l|Leventhal. MacCallum. fc Stanj Il978t ) . Sin ce then many 



at the centre of t h e Milky Way with a FWHM of 6-8° (e.g 
Jean et al.l l2003l : iKnodlseder et al. I l2005l : IChurazov et all 



ballo o n flights and several space missions (e .g. iGehrels et al 



199ll: iMahonev. Ling, fc Wheatoiil 1 1994 IXeegarden et al 



19961 : iPurcell et al.l 199?! ) have measured the spatial distri- 
bution and spectral properties of the line. 

Since 2003, INTEGRAL data on the 511 keV line emis- 
sion became available. These data show that there is a 
strong, 10~^phot s~^ cm~^, source of 511 keV photons 



while the 511 keV flux coming from the Galactic 
plane is much less constrained by the INTE GRAL data (e.g. 
iTeegarden et al.ll2005l : iBouchet et al.l 120081 ). No significant 
variability in the narrow 511 keV line ha s been found in the 
SPI data (|Tsvgankov fc Churazovll2010l ). The spectrum of 
the annihilation radiation from the Galactic Centre region 
observed by SPI can be reasonably well described as a com- 
bination of a narrow line (~2-3 keV Gaussian at 511 keV) 
and a three-photon continuum. The flux ratio between the 
three-photon continuum and the narrow line suggests that 
the majorit y of positrons form positronium prior to annihi- 
lation (e.g. Ichurazov et~al.,, 2005. : , Jean et al.. 20061 ). 

In spite of numerous observations, the origin of the 511 
keV emission from the GC is not established. The problem 
can be summarised as follows: 

• The morphology of the observed 511 keV line emission 
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does not fit the disk-hke spatial distribution of such obvious 
sources of positrons as i) massive stars, which produce /J"*" 
unstable nuclei like ^®AI or ii) tt^'s generated by the inter- 
action of cosmic rays with the interstellar medium (ISM). 
This favours scenarios with a more prominent central excess 
("bulge dominated"), including the old stellar population 
(e.g. supernovae la), activity of Sgr A* and the annihilation 
of dark matter particles. 

• The spectrum of the annihilation emission (narrow line 
width and the flux ratio of the line and three-photon con- 
tinuum) suggests that positrons are annihilating in a warm, 
~ 10^ K, and slight l y ionized, ^ a f ew %, medium (e.g. 
IChurazov et al.ll2005l : fjean et al.ll2006l '). Such an ISM phase 
does exist in the Galaxy, but its spatial distribution is 
strongly concentrated towards the plane, posing the problem 
of "finding" the right ISM phase by a positron. 

Recently an atte mpt to build a self-consistent pii 



ture was presented by iHigdon. Lingenfelter. fc Rothschild! 
(|2009h : lLingenfeher. Higdon. fc Rothschild! (|2009h . who as- 
sumed that the spatial propagation of positrons, pro- 
duced via /3+ decay of ^^Ni, ^^Ti and ^"^Al, is governed 
by a diffusion process with the effective diffusion coeffi- 
cient different in the bulge and the disk of the Galaxy. 
During diffusion the positrons enter the HII and HI en- 
velopes of m olecular clouds, in particular t hose forming a 
"Tihed Disk" jFerriere. Gillard. fc Jeanl!2007!') within 1.5 kpc 



of the GC. The model of lHigdon. Lingenfelter. fc Rothschild! 
(!2009l ). with a reasonable set of assumptions, can explain the 
basic properties of the annihilation radiation. 

Here we take an alternative route and consider a pos- 
sibility that a significant fraction of positrons are born in 
the hot ISM, which is eventually able to cool via radiative 
losses and perhaps via adiabatic expansion. We show below 
that with these assumptions the spectral properties of the 
511 keV line can be easily explained. However, to make our 
picture self-consistent one needs to model the thermal state 
of the ISM, which is beyond the scope of the present paper. 

This paper is based on the INTEGRAL data accumu- 
lated over ~ 6 years of observations and aims at placing 
tighter constraints on the spectral and spatial properties of 
the annihilation emission. 



2 DATA AND BACKGROUND HANDLING 

SPI is a co ded mask germanium spectrometer on board IN- 
TEGRAL (! Winkler et al.!!2003l ). launched in October 2002 
aboard a PROTON rocket. The instrument consists of 19 
individual Ge detectors, has a field of view of ~30° (at zero 
response), an effective area at 511 keV of ~ 70 c m^ and en- 
ergy resolution of ~2 keV jVedrenne et al.l!2003! ). The good 
energy resolution makes SPI an appropriate instrument for 
studying the spectrum of e+e~ annihilation emission. 

For our analysis we use all data available to us by mid- 
2009, including public data, some proprietary data (in par- 
ticular, proposals 0420073, 0520071 and parts of 0620059). 
Prior to actual data analysis, all individual observations 
were screened for periods of very high particle background. 
We use the SPI anticoincidence (ACS) shield rate as a main 



indicator of high background. Several additional observa- 
tions were also omitted from the anal ysis, e.g. those take n 
shortly after SPI annealing procedures (|Roaues et al.ll200^ '). 
For our analysis we used a combination of single and pulse- 
shape-discriminator (PSD) events fsee !Roaues et al] !2003! . 
for details) and treated them in the same way. 



2.1 Energy gain and resolution 

For each detector, a linear relation between the energy and 
the channel number was assumed and calibrated (separately 
for each orbit), using the observed energies of lines at ~ 
198, 438, 584, 882, 1764, 1 779, 2223 and 2754 keV (see 
!Weidenspointner et al.!l2003l . for a comprehensive list of SPI 
background lines). With this calibration the RMS deviation 
of the background 511 keV line energy (revolution based) is 
0.0066 keV while the mean energy of the line is 510.926 keV. 
There is thus a small deviation of the line mean energy from 
the electron rest energy (510.999 keV): AS = 0.07 keV; 
it can be attributed to the simplified linear energy /channel 
relation. This deviation is comparable to the statistical un- 
certainty on the line energy (see ^ and no attempt was 
made to correct for this effect. 

Since the background 511 keV line (produced by the 
positrons annihilating in the body of the detector) is kine- 
matically broadened, we used the two bracketing lines (at 
438 and 584 keV) to calculate the resolution at 511 keV as 
FWHMsii = 0.5 X (FWHM438 + FWHM584). The resuh- 
ing value is FWHM511 = 2.175 keV when averaged over all 
observations in the vicinity of the Galactic Centre. 



2.2 Background modeling 

In the background mo deling we followed the scheme used in 
!Churazov et al.! (!2005! ). Namely, the background count rate 
B{i,E) in detector i at energy E is assumed to be propor- 
tional to the detector saturated (i.e. above 8 MeV) event 
rate Rsat and time t: 



B(i, E) = C{i, E) + a{i, E) x i?,at + /?(i, E) x t, 



(1) 



where i = 1, 19 is the detector number. The coefficients 
a(i,E) and /3{i,E) of this linear relation were determined 
separately for each detector/energy channel using all the 
available SPI data, while the constant C{i, E) was estimated 
using only the data away from the Galactic plane and away 
from bright sources. The whole data set was divided into 
14 time intervals (defined by the annealing periods and the 
dates of individual detector failure^) and the coefficients 
a{i,E), l3(i,E) and C{i,E) were determined separately for 
each interval. This procedure, although not providing a per- 
fect description of the background, works reasonably well at 
all energies and, given the small number of free parameters, 
is very robust. For instance, the relative RMS deviation of 
the 600-1000 keV fiux (averaged over one revolution) from 
the count rate predicted by this model is ~ 1.3 10"'^. 

Along with our reference background model, we made 



1 Tliree out of the 19 SPI detectors liad failed by tiie end of 2009. 
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Figure 1. Exposure map (in Galactic coordinates) showing all the observations used in the present analysis. Off-Galactic- plane obser- 
vations were used for determining the instrument background. For each observation, the effective observing time is added to the pixel in 
the map containing the pointing direction of the telescope's axis. The total exposure is about 70 Ms. Pixels are 2 degrees on a side. 



several experiments by dividing the full data set into differ- 
ent time intervals and modifying the selection criteria used 
for cleaning the data from background flares. In these exper- 
iments we generated several "alternative" background mod- 
els. Below we use these models to test the sensitivity of our 
results to the details of the background modeling. 



3 IMAGING 

Imaging with coded mask telescopes often requires addi- 
tional assumptions to be made on the sky surface brightness 
distribution, e.g. sparsity of compact sources or regular be- 
haviour of diffuse emission. This is especially important for 
the SPI telescope, which has only n — 19 independent de- 
tectors/pixels. This makes image reconstruction in a single 
observation problematic, given that a typical image recon- 
struction procedure relies on the smallness of the parameter 
1/n with respect to unity. We have chosen not to make a 
priori assumptions about the spatial distribution of sources 
and restricted ourselves to three different types of "imag- 
ing" : (i) light bucket imaging, (ii) simple parametric models 
and (iii) linear decomposition of the sky surface brightness 
into few simple templates. 

3.1 Light bucket mapping 

In this simple "imaging" procedure, for each observation, 
the background subtracted counts in a given energy band 
are summed over all detectors and the resulting signal is 
added to the pixel of the map corresponding to the telescope 
pointing direction. Similarly the exposure map is formed by 
adding the exposure times (see Fig. [ij. The ratio of the 



counts and the exposure map gives the count rate image, 
which is crudely converted into units of phot/s/cm^ using 
the effective area for an on-axis point source. Such maps will 
provide a correct point source flux if the telescope is pointed 
directly towards the source. But given the large size of the 
SPI FoV (~30° diameter at zero response) and sparse cover- 
age of the sky, 'light bucket mapping' produces rather crude 
maps, while the procedure itself is very simple and robust. 
We use hght bucket mapping to get a global surface bright- 
ness map or to follow the variations of the surface brightness 
along and across the Galactic plane with poor angular res- 
olution. Note that in this mode SPI is used as a coUimated 
instrument and its mask imaging capabilities are not ex- 
ploited. The net flux (difference between the observed count 
rate and predicted background model) crucially depends on 
the quality of the background model. 

3.1.1 Scans along the Galactic plane 

In Fig. [2] we show slices of the Galaxy along the Galac- 
tic plane in three energy bands. Each slice has a width of 
16° perpendicular to the Galactic plane and is centred at 
6 = 0°. The step along the Galactic plane is 2° for the two 
upper panels and 4° for the bottom panel. 

In the upper panel (50-100 keV), three prominent peaks 
correspond (from left to right) to Cyg X-1, Galactic Centre 
and the Crab Nebula. The width ~30° and the complex 
structure of the peaks near the maximum (clearly visible for 
Cyg X-1) are due to the presence of the mask and partly 
due to the intrinsic variability of the sources. The few re- 
gions with negative fluxes correspond to observations where 
the actual background count rate is higher than the model- 
predicted background described in ^2.2\ Note that the model 
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Figure 2. Scans (light bucket imaging) along the Galactic plane in 3 energy bands. Fluxes are averaged over 2° bins over I and ±8° over 
b. For the 1804-1813 keV band, the width of the bins along i is 4° . By construction all points are statistically independent (except 
through the background subtraction). Vertical grey bars show the positions of the series of scans perpendicular to the Galactic plane 
discussed in i|3.1.2l The positions of the Crab Nebula, Cyg X-1 and Centaurus A are marked with the vertical dotted lines. 



background was calculated using a set of "blank fields", 
whose definition is problematic given the large size of the 
SPI FoV. The sets of blank fields are different in each of 
the 14 time intervals used in the background modelling (see 
tj2.2[) . For this reason, it is possible that the actual flux from 



the "blank fields" is higher than that from some patches of 
the Galactic place free of strong sources. 

The morphology of the slice in the 508-514 keV bantfl, 

^ See iTeeearden et al" I 1I2OO5I ') for an earlier version of a similar 
plot. 
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containing the 511 keV line (Fig. [2l middle panel), is 
markedly different. With the 2° bins along the Galactic 
plane, the only prominent feature is the peak at the Galactic 
Centre, with the extent along I roughly similar to the extent 
of the peak near Cyg X-1 at low energies. The peak flux is 
~ 10~^ phot s~^ cm~^. As already emphasized, this value 
is not a precise measure of the true 511 keV flux from the 
Galactic Centre region, but it should be accurate to a factor 
of better than 2 if the source size does not exceed several 
degrees. 

No strong asymmetry of the positive and negative 
longitude wings of the central pe ak is apparent in the 
light bucket profiles (see, however, IWeidenspointner et al.l 



12008 , who report an ex cess flux for negative longitudes 



Bouchet et al.l [20081 . who do not find strong evidence 



and 

for asymmetry). An excess at negative longitud es I ~ - 
25° is visible in the light bucket profile in Tccga rden et al.l 
(|2005h (their Fig. 2). Renormalization of their fluxes to 
the units used in Fig. [5] suggests an excess at the level 
~ 1.5 10^^ phot s^^ cm~^. In our analysis the flux at / ~ 
-25° is essentially consistent with zero. If anything, the pro- 
files shown in Fig. [2] suggest an excess in the 508-514 keV 
flux at positive longitudes / ~ 25° . However, our experi- 
ments with "alternative" background models mentioned in 
tj2.2l showed that a weak spurious asymmetry (corresponding 
to a flux excess/deficit of ~ 10"'* phot s~^ cm~^ at |/| ~20- 
25° ) can appear at either negative or positive longitudes. 
We therefore conclude that, with the present knowledge of 
the SPl background, light bucket profiles do not provide 
compelling evidence for asymmetry in the 508-514 keV flux 
along the Galactic plane. 

Finally in the bottom panel of Fig. (2] the slice of the 
Galactic plane in the 1804-1813 keV band, containing the 
1.8 MeV line of ^''Al, is shown. Unlike the 511 keV emission, 
there is a broad peak (much broader than the SPI response 
to a point source), centred at the GC. This distribution is 
quali tatively consistent with the Co mptel (jPliischke et all 
bOOll ) and earlier INTEGRAL resuhs (|Wang et al.ll2009l 'l. 



3.1.2 Scans perpendicular to the Galactic plane 

A similar approach of light bucket mapping can be applied 
to the special scans of the Galactic plane along b performed 
in 2007, 2008 and 2009. The positions of the scans in / are 
denoted with the gray vertical lines in Fig. (2] Three out of 
four scans were made as a sequence of pointings starting and 
ending ~ 30° away from the plane on either side of it. The 
pointing direction was changing in ~2° steps. The remain- 
ing "scan" at I — had been done earlier as a set of point- 
ings directly at the plane and ~25° above and below the 
plane. The advantage of making "fast" (taking 10-15 hours) 
scans over a selected region of the sky is that the quality of 
the background modeling can be verified/improved if typical 
variations of the detector background occur on longer time 
scales. 

The scans in the 508-514 keV band in the direction 
perpendicular to the Galactic plane at i ~ 0° , 22° , -45° are 
shown in Fig. [3] The signal is clearly visible at I ~ 22° and is 
close to zero for I ~ -45° . To guide the eye, we have plotted 



with the black curves the expected light bucket profiles for 
the case where the spatial distribution of the flux in the 508- 
514 keV band follows the "Bulge" model described later in 
>]3.5I The expected profiles were crudely estimated assuming 
that the axis of the instrument was moving along b at fixed 
I, while keeping the position angle fixed. This causes visible 
asymmetry in the peak of the I ~ 0° profile, reflecting the 
particular orientation of the SPl coding mask relative to the 
source. In reality, the observed profiles are combinations of 
observations with varying position angle and varying I and 
hence they may differ from the simulated profiles in subtle 
details. Broadly, the light bucket imaging shows the overall 
consistency with a bulge-dominated distribution of the 508- 
514 keV fiux and indicates that the disk emission is weak at 
/ ~ -45° . 



3.2 Simple parametric models 

We now use a simple function - two-dimensional Gaussian 
- to model the spatial distribution of the annihilation line 
emission near the Galactic Centre. The data collected when 
the INTEGRAL pointing direction was within 30° of the 
GC were used. Two models were considered: 



L^ln2.b±l 



Gil, b) = Fi X e I 



GD{l,b) =FiXe I "1^ 



+ F2X e 



(2) 



where G{1, b) is the surface brightness distribution as a func- 
tion of Galactic longitude and latitude, F is the flux and Wi, 
Wb are the full widths at half maximum along I and b, re- 
spectively. The second model has an additional component 
that is aimed to account for emission elongated over the 
Galactic plane. Since only the data within 30° of the GC 
were used (30° corresponds to the deviation of the SPl axis 
from the GC), the "infinite" extent of this component over 
/ means that the surface brightness of this component does 
not decrease much at a distance of ~45° from the GC. 

To verify the sensitivity of the results to a particular 
functional form, we also used an exponential law instead of 
the Gaussian: 



Eil,b) 
ED{l,b) 



Fi xe 
Fi X e 



(3) 



For a given pair of Wi and W h, the model is conv olved 
with the simulated SPl response jSturner et al.ll2003l ) and 
compared with the count rate in the 508-514 keV band in in- 
dividual detectors during individual observations. The best- 
fitting value of Fi for the one-component model or the pair 
of values Fi and F2 for the two-component model is calcu- 
lated. Here and in all the subsequent analysis, we use simple 
statistics. The problems of using the criterion for low 
photon c ounting statistics are circumvented following the 
recipe oflCh urazov et al.l l)l996h . Namely, the standard de- 
viation associated with the count rate in a given spectral, 
spatial or time bin is evaluated using the mean count rate 
averaged over a large number of "nearby" similar bins. 

Typical values of x'^ per degree of freedom for our mod- 
els are ~ 1.01. Given the large number of degrees of free- 
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Figure 3. Scans (light bucket imaging) perpendicular the Galactic plane at i ~ 0° , 22° , -45° in the 508—514 keV energy band. The 
solid and dashed curves approximately show the expected profiles at Z ~ 0° and 20° if the spatial distribution of the flux in the 508—514 
keV band is adequately described by the "bulge" component described in i|3.5l 



dom (~ 250000 = number of detectors times the number of 
observations) and the very low signal-to-noise ratio of the 
annihilation signal in individual observations, this value is 
not a useful indicator of the "absolute" quality of the model. 
Instead, the change of can be used to compare different 
models or place constraints on the model parameters. 

The resulting contours of for both models (as a func- 
tion of the parameters Wi and Wb) are plotted in Fig. U 
The dashed lines correspond to the one-component model 
(pure 2D Gaussian), while the solid lines correspond to the 
two-component model GD{l,h) (see eq. [2}. The black dots 
mark the positions of 'x^ minima. Contours are spaced by 
Ax^ = 3. The best-fitting values of the widths are 10° x 
6° (if the one-component model is used) and ~8° x 6° (for 
the two-component modeljf]- 

Clearly, for the pure 2D Gaussian model the data sug- 
gest a significant flattening of the distribution towards the 
plane. If an additional component, extended along the plane, 
is included (the two-component model in eq. then the 
best-fitting central Gaussian is much more symmetric in I 
and h. The improvement in x'^ for the two-component model 
compared to the one-component model is ~ 40. The position 



^ Note that we used a grid over FWHM values with a 1° step and 
a 7° X 6° Gaussian gives almost the same ^ the best-fitting 
8° X 6° Gaussian. 



of the minimum for the two-component model does not de- 
pend much on the width of the extended component over 6. 

We tried for the second component Wd = 2° , Wd ~ 6° and 

I IM ] 

Wd ~ 10° as well as an exponential shape e 1^ ™d I instead 
of the Gaussian, and got essentially the same best-fitting 
parameters (8° x 6° ) for the central Gaussian. 

We emphasize here that the presence of the second com- 
ponent does not necessarily imply that the Galactic disk 
is "detected", but rather that a single symmetric Gaus- 
sian/exponential component is not a perfect description of 
the data. 

f I'l] f IM ] 

An exponential shape e 1 J I with Wi ~3° and 

Wb ~2° can also be used to describe the central component 
(see Table [1} instead of a Gaussian. In fact, the largest im- 
provement in the is achieved when both the central and 
disk components are described as exponential functions (Ta- 
ble [l|. While the difference in the x^ between various mod- 
els (e.g. one-component exponential versus one-component 
Gaussian) is formally statistically significant, it corresponds 
to less than 1% change in the Ax'^ value relative to the "null 
hypothesis" of a zero flux in the 508-514 keV band over the 
entire data set. At this level, a particular value of Ax^ for a 
given model might be sensitive to the subtle features of the 
data set. Our experiments with different background models 
have indeed shown that the ranking of models in terms of 
X^ can slightly vary. For example, the ranking of the Gaus- 
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Figure 4. Contours of ^ function of the FWHM along 

I and b when the spatial distribution of the 508-514 keV flux 
is described by a 2D Gaussian. The dashed lines correspond to 
the pure 2D Gaussian model {G{l,b) in cq. [2)l, while the solid 
lines correspond to the model Gd{1, b) in eq. [2] which includes an 
additional extended component. Black dots mark the positions of 
the minima. Contours arc spaced by Ax^ = 3. The best-fitting 
values of the width arc 10° X 6° and 7-8° X 6° for the first and 
second model, respectively. 



sian+Gaussian and Gaussian+exponential models in Table 
[T]can change depending on the particular background model. 

Despite the difference in the functional forms, the 
flux in the central component does not depend strongly 
on the model used; it varies between 8.4 10~* and 
9.3 10~* phot s~^ cm~^ for all two-components models. 



3.3 Position of the centroid 

Using the one-component exponential model (the first model 
in eq. [3]) with Wi = 3° and Wb — 2° , we tried to vary the 
centroid of the distribution over a few degrees around the 
GC and calculated the changes in the x^- The best-fitting 
position is at {I, b) = (—0.1°, —0.2°) and the minimal dif- 
fers from the value at {I, b) = (0, 0) by 2.3. While this means 
that formally the statistics of accumulated data is already 
sufficient to measure sub-degree shifts of the centroid, we 
believe that systematic uncertainties and the obvious crude- 
ness of the spatial model preclude any firm conclusion. We 
can conservatively conclude from this exercise that the po- 
sition of the centroid is consistent within ~0.2° with the 
position of the dynamic centre of the Milky Way. 



Figure 5. x^ ^ function of the tilt of a two-dimensional Gaus- 
sian (solid black line). The tilt is measured relative to the Galac- 
tic plane. Positive values of the tilt correspond to the major axis 
pointing towards positive b at negative I. For each value of the tilt 
the best fitting values of the FWHM in two dimensions were iden- 
tified on a crude FWHM grid with steps of ~1° . The ratio of the 
FWHMs (Wi/Wii) is shown in the bottom panel. The minimum 
of the is achieved at 61 = 12° and Wi = 14° , VKj, = 4° . The 
large variations of the Wi /Wb ratio are owing to the crudeness of 
the grid and a very shallow minimum in the function X^i^l/^b) 
for a given 9. As \9\ increases the Gaussian converges to a sym- 
metric structure {Wi/Wb ~ 1). For comparison, the red dotted 
line shows the x^ as a function of the tilt if Wi and Wb are fixed 
at their best-fitting values of 10° and 6° , respectively, for 6 = 0. 



3.4 Tilted Gaussian 



iHigdon. Lingenfelter. fc Rothschildl (|2009l ') have suggested 
that a significant fraction of positrons is annihilating in a 
"Tilted Disk" of neutral gas inside the central 3 kpc re- 
gion of the Galaxy. This "Tilted Disk" is one of the compo- 
nents of the iFerriere? Gillard. fc jeani (|2007D model of the 
interstellar gas distribution in the innermost part of the 
Milky W ay, based on earlier res ults of lLiszt fc BurtonI (|198G| ) 
fs ee alsolLiszt fc Burtonlll996l'l. T he Tilted Disk (see Fig. 4 
in iFerriere. Gillard. fc Jean in projection to the sky 



plane has an apparent size of ~18° by ~5° and is tilted by 
~ 30° with respect to the Galactic plane. Since the extent 
of the Disk in I and b resembles the dimensions of th e 511 
keV source, iHigdon. Lingenfelter. fc Rothschildl (|2009l ) sug- 
gest that the positrons reach the disk and annihilate there. 

To test further the Tilted Disk model we fit the data 
with a two-dimensional Gaussian allowing for rotation of 
the major axis. The results are shown in Fig. [5] The best- 
fitting rotation angle and the widths of the Gaussian in two 
directions are ~12° and Wi = 14° , Wb = 4° , respectively. 
The improvement in compared to = 0° (and Wi = 10° , 
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Figure 6. The model templates used for fitting the data in the 
508-514 keV band (from top to bottom): Disk, Disk with a Hole, 
Bulge. The projected surface brightness distributions shown arc 
base d on the 3D models of stellar density distribution ad opted 
from lL6pez-Corredoira. Cabrera-Lavers. fc Gerhar i Hooi). 

Wb = 6° ) is Ax'^ ~ 7. As discussed already, the values of 
Ax'^ of the order of a few are formally statistically significant 
if the errors are due to counting Poisson noise only. The total 
change of the when a two-dimensional Gaussian is added 
to the model is ~ 8000. Comparing these numbers, it is clear 
that even modest systematics in the data could affect precise 
derivation of the source spatial characteristics. Taking the 
results shown in Fig. [5]at face value, a marginally significant 
improvement in the fit is possible if the annihilation region 
is tilted by ~ 12° with respect to the Galactic plane. 

As \6\ increases beyond 25° , the Gaussian con- 
verges to an almost symmetric structure. The tilt of ~ 
30° corresponds to Ax^ ~ 23 for an almost symmetric Gaus- 
sian with Wi = 9° and Wb = 7° . We therefore conclude that 
a vers ion of a the Tilted Disk of lFerriere. Gillard. fc JeanI 
(|2007h is not particularly favored by SPI data compared to 
a structure with zero tilt. However, moderate values of the 
title <20° are allowed by the SPI data. 

3.5 Decomposition using plausible templates 

The simple hght bucket imaging and parametric fitting of 
the 508-514 keV flux done in i|3.1l and l3.2l suggest that a rea- 
sonable description of the data in the central radian could 
be achieved with an almost symmetric Gaussian at the cen- 
tre and a more extended component along the plane. The 
stellar bulge and stellar disk are the most obvious structural 
components of the Galaxy that qualitatively fit this descrip- 
tion. In fact, there are plenty of disk-like structures (cold 
gas, massive stars, cosmic ray induced gamma-ray flux, etc) 
and very few components that show a prominent peak at 
the GC. One of the known centrally-peaked distributions 
is that of the NIR hght, which is a tracer of the old stel- 
lar population of the Galaxy. For this reason, we decided 
to restrict our analysis to three simple templates, corre- 
sponding to stellar components of the Milky Way, adopted 



from iLopez-Corredoira. Cabrera-Lavers. fc Gerhard! (|2005l ) 
and shown in Fig. [G] 

• Disk (eq. 2-5 in Lopez-Corredoira et al.) 

• Disk with a Hole (eq. 6-7 in Lopez-Corredoira et al.) 

• Bulge (see §3 in Lopez-Corredoira et al.). 

This is by no means a comprehen s ive li st of possible tem- 
plates (see e.g. iKnodlseder et akl |200^ for a systematic 
analysis of various templates). However, the hope is that 
these spatial templates will capture the most basic proper- 
ties of the 508-514 keV flux distribution, even if the physical 
motivation is questionable or the details of shape are not 
correct. 

We also used slightly modified (truncated) versions of 
these templates; for each template we set to zero the surface 
brightness in all regions where it is smaller than 10% of 
the peak value. The motivation behind this modification is 
an attempt to have a template that has a bulge or disk- 
type morphology at high surface brightness but lacks the 
extended low surface brightness regions, which might make 
a significant contribution to the total flux. 

For a given set of templates, the surface brightness for 
each template is con volved with the simulated SPI response 
(|Sturner et al.|[2003l ') yielding an expected count rate in the 
508-514 keV band in individual detectors during individ- 
ual observations. The best-fitting normalizations of the tem- 
plates are then calculated in order to minimize the devi- 
ation between the raw data and the model. 

The x^ values for the various spatial models are given 
in Table [T] The fluxes quoted in Table [T] are the integrated 
model fluxes (with best-fitting normalization) over a square 
80° X 80° region around the GC. The choice of the region 
for the calculation of the integrated model fiux is rather ar- 
bitrary. The 80° X 80° region is not uniformly covered by 
observations, and areas with the smallest errors (largest ex- 
posures) dominate in the determination of the best-fitting 
normalization. On the other hand, underexposed (or even 
not observed at all) areas can still provide a significant con- 
tribution to the total integrated flux. An example of a clear 
overestimation of the flux due to this effect is seen for the 
pure "Disk" models in Table [T] These models have poor 
, but predict a large integrated flux, since their normal- 
ization is largely set by the innermost bright region of the 
Galaxy. This becomes an especially severe problem when 
dealing with fluxes integrated over very large areas (e.g. the 
flux from the entire Galactic disk). For this reason, we quote 
fluxes integrated over the 80° x 80° region rather than from 
the whole Milky Way. 

Table [1] strongly suggests that the spatial distribution 
of the 508-514 keV flux in the inner region of the Galaxy 
is more extended along the Galactic plane than perpendic- 
ular to it. This is clear from the parameters of the single- 
component parametric models (see also Fig. [4]) and from 
the comparison of the x'^ values for the one-component and 
two-component (e.g. two Gaussians) models. 

Among the one-component models, the exponential 
model has a slightly better x^ than the Gaussian model or 
the Bulge model. For the two-component models, the expo- 
nential Bulge -I- Gaussian Disk has the smallest x^ among 
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Table 1. Fitting the 508-514 keV surface brightness with templates and parametric models. Quoted Ax^ values were calculated relative 
to the "null hypothesis" of zero flux in the 508-514 keV band over the entire data set. The data within 30° from the GC were used. 
Quoted fluxes correspond to the integrated model fluxes (with best-fitting normalization) over a 80° X 80° square region around the 
GC. 



Spatial template — Ax^ Fi F2 

Component 1 Component 2 phot s^^ cm^^ phot s^^ cm~^ 



Templates described in i|3.5l with the surface brightness truncated at 10 ^ of majcimum 



Disk 




6837.7 


(2.20 ±0.03) 10" 


3 






Disk/Hole 




6660.3 


(2.22 ±0.03) 10- 


3 






Bulge 




8252.4 


(1.12 ±0.01) 10- 


3 






Bulge 




Disk 8252.5 


(1.11 ±0.03) 10- 


3 


(2.4 ±6.4) 10" 


'S 


Bulge 




Disk/Hole 8252.5 


(1.11 ±0.03) 10- 


3 


(2.0 ±6.2) 10- 


-5 




Templates » 


ith the surface brightness truncated at 10 ^ of maximum 






Disk 




6964.0 


(1.83 ±0.02) 10- 


-3 






Disk/Hole 




6758.9 


(1.85 ±0.02) 10- 


3 






Bulge 




8229.1 


(0.97 ±0.01) 10- 


3 






Bulge* 




Disk 8252.9 


(0.87 ±0.02) 10- 


3 


(2.4 ±0.5) 10" 


-4 


Bulge 




Disk/Hole 8253.0 


(0.87 ±0.02) 10- 


3 


(2.3 ±0.5) 10- 


-4 


G(l,b) and G]j(l,b) models from eq. [2]and[3] 


Gaussian Wi =10 


' , Wb =6° 


8223.4 


(0.96 ±0.01) 10- 


3 






Gaussian Wi =8° 


, Wb =6° 


Gaussian, Wd=2° 8260.3 


(0.84 ±0.02) 10- 


3 


(2.9 ±0.5) 10- 


-4 


Gaussian W; =8° 


, Wb =6° 


Gaussian, Wd=6° 8262.8 


(0.84 ± 0.02) 10- 


3 


(3.1 ±0.5) 10" 


-4 


Gaussian W; =8° 


, Wb =6° 


Gaussian, Wd=10° 8263.7 


(0.84 ±0.02) 10- 


3 


(3.3 ±0.5) 10" 


-4 


Gaussian Wi =8° 


, Wb =6° 


Exponential, 1^0=2° 8262.6 


(0.84 ±0.02) 10- 


3 


(3.1 ±0.5) 10- 


-4 


Gaussian Wi =8° 


, Wb =6° 


Exponential, 1^0=6° 8262.6 


(0.84 ±0.02) 10- 


3 


(3.7 ±0.6) 10- 


-4 


Exponential Wi = 


3° , Wb =2° 


8265.4 


(1.01 ±0.01) 10- 


3 






Exponential Wi = 


3° , Wb =2° 


Gaussian, 1^0=2° 8284.1 


(0.93 ±0.02) 10- 


3 


(2.1 ±0.5) ID- 


-4 



* - The Bulge component of the Disk+Bulgc model is used in the spectral analysis in J5] 



the set of models considered here. It is clear that i) the true 
distribution of the 508-514 keV flux is likely more compli- 
cated than any of the models used and ii) the relatively 
modest (although statistically significant) changes in the 
among different models (containing the Bulge component) 
suggest that the basic properties of the distribution are cap- 
tured by our models. 

The uncertainty in the choice of the spatial model 
directly translates into the uncertainty of the 508-514 
keV flux. In particular, the total flux of the Bulge com- 
ponent varies among all two-component models (both 
parametric and based on templates) from 0.84 10"'' to 
1.11 10""^ phot s~^ cm~^. At the same time, the total Bulge 
+ Disk flux (integrated over the 80° x 80° region) varies 
from 1.10 10"^ to 1.14 10"^ phot s"^ cm"^ 

Since some of the spatial templates used in Table [1] 
are projections of the 3D models (based on stellar distri- 
bution), it is easy to recalculate the observed flux F into 
the total luminosity L of the corresponding component as 
L = A-kD^ X F. Here D depends on the 3D distribution of 
the volume emissivity and on the region used to calculate 
the flux F. The values of D are given in Table (2] Using 
the effective distances from Table [2] for the templates trun- 
cated at 10% of the maximal surface brightness, the total 
luminosity of the Bulge component is ~ 5.0 10*^ phot s~^, 
while the total Disk luminosity is ~ 2.3 10*^ phot s~^. As 



Table 2. EfEcctivc distance D needed to recalculate the observed 
flux from the 80° X 80° region around the GC into the lumi- 
nosity for different spatial components. Two columns correspond 
to different cutoffs in the surface brightness distribution of each 
component (relative to the maximal surface brightness of this 
component). 



Spatial template 


Effective distance, kpc 




cutoff=10-5 


cutoff=0.1 


Disk 


7.74 


8.98 


Disk/Hole 


7.79 


8.99 


Bulge 


5.89 


6.93 



noted above, the resulting luminosities rely on the model 3D 
distribution and for this reason they are model dependent. 



4 ^'■Al DECAY 

^®A1 is one of the obvious and ac countable sources of the 
positrons in the Galaxy (see e.g. iDiehl et al.ll2003 ). This 
radioactive isotope with the half-life time of 7.17 10^ years 
is produced by massive stars, which primarily occupy the 
disk of the Galaxy. The positrons are produced in 81.7% 
of decays, with the mean energy of 543 keV; in 99.8% a 
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1.809 MeV photon is emittec(f|. Assuming that the fraction 
of positrons forming positronium is /ps (see ij5}, the fiuxes 
in the 511 keV and 1.8 MeV lines per 1 positron produced 
via ^^Al decay are: 

F511 = 0.817 [(1 - /p.) + i/p.] X 2 = 1.63 [1 - |/p.] , 
Fi.s = 0.998. 

Here we assume that i) the fraction of para-positronium is 
1/4 and ii) all positrons annihilating without formation of 
positronium (fraction 1 — /ps) produce a narrow 511 keV 
line. The factor 2 in the above expression accounts for the 
2 photons produced in two-photon annihilation. Therefore, 
F511 = 0.409, 0.471 and 1.637 x Fi.s, for /p^ =1, 0.95 and 0, 
respectively (see 33 for details). 

It is interesting to compare this relation with the fluxes 
obtained from the light-bucket images in the 508-514 keV 
and 1805-1813 keV energy bands, which should contain 
most of the 511 keV and 1809 keV line fluxes, unless the lines 
are strongly broadened. Shown in Fig.[7]is a light-bucket lon- 
gitude scan of the Galactic plane over the region I = ±60° . 
The blue thick crosses show the flux in the 508-514 keV 
band, while the thin red crosses show the observed flux in 
the 1805-1813 keV band scaled by a factor of 0.409. These 
red crosses correspond to the expected 511 keV line flux 
arising from '^^Al decay under the assumption that 100% 
of positrons annihilate through the formation of positron- 
ium. It will be observed later (^ that this assumption is 
supported by spectral data. Clearly, the high 511 keV line 
flux from the GC region cannot be explained by ^^Al de- 
cay (unless the positrons produced in the disk are someh ow 
transported to the central region; see e.g. IPrantzodlioO^ V 

One can use the same set of templates as in tj3.5l and 
repeat the analysis for a set of energies to obtain the spec- 
trum associated with each spatial component. We did this 
using the Disk-|-Bulge model, which was applied to all ob- 
servations with the INTEGRAL pointing direction within 
30° from GC. Since we do not explicitly account for com- 
pact sources in this procedure, the low energy parts of the 
spectra obtained may not be reliable. However, near the 511 
keV and 1.8 MeV lines the contribution of individual sources 
should not be crucial. Nevertheless, this simplifying assump- 
tion should be kept in mind when interpreting the resulting 
spectra, shown in Fig. [S] 

The upper row in Fig. [S] shows the spectrum of the 
Bulge component in the vicinity of the 511 keV and 1.8 MeV 
lines. The red line in the top-left panel shows the best-fitting 
model (Gaussian at 511 keV and ortho-positronium contin- 
uum; see IJ5] below for details). The fiux in the 511 keV line 
in this model is F511 = (0.84 ± 0.03) 10"^ phot s~^ cm"^. 
Note that this flux was obtained from the spectral analysis 
and it differs slightly from the "Bulge" flux obtained from 
the analysis of the count rate in the 508-514 keV band, 
quoted in Table [T] The value F511 was then recalculated 
into the expected flux in the 1.8 MeV line using the relation 
Fi,s = F511/O.409, corresponding to the assumption that all 
positrons are produced by ^^Al decay and the fraction of 
annihilations via positronium is 100% (see eq. [4]and Sj5]). A 

* http://www.nndc.bnl.gov/mird/. 
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Figure 8. Spectra of the Bulge (top row) and Disk (bottom 
row) near 511 kcV and 1809 keV, obtained using the decompo- 
sition of the data with the two-component Disk-|-Bulge model. 
The red line in the top-left panel shows the best-fitting model 
(Gaussian at 511 keV and ortho-positronium continuum; see ijS] 
below for details). The flux in the 511 keV line in this model 
is -F511 = 0.84 10"'' phot s~^ cm~^. The blue curve in the top- 
right panel shows a Gaussian line at 1.809 MeV. The flux in the 
Gaussian was calculated as Fi,s = -F511/O.409, corresponding to 
the assumption that all positrons are produced by ^®A1 decay 
and the fraction of annihilations via positronium is /pg = 1.0 
(see eq. |4] and ijS) . The bottom row shows the spectrum of the 
Disk component. In contrast to the Bulge spectrum, the 1.8 MeV 
line is very prominent in the spectrum (bottom-right panel). The 
red line shows the best-fitting Gaussian at 1.809 MeV, with flux 
Fi.g = 4.1 10^''phot s~^ cm~^. The blue line in the bottom-left 
panel now shows a 511 keV line with flux F511 = 0.409 X Fi.g. 



line with the resulting flux Fi g is plotted in the top-right 
panel as a Gaussian line centered at 1.809 MeV. Clearly 
such a strong line at 1.8 MeV is in stark contrast to the 
data, which show no evidence for 1.8 MeV emission in the 
Bulge spectrum. 

The bottom row in Fig. |8] shows the Disk component 
spectrum in the vicinity of the 511 keV and 1.8 MeV lines. 
In contrast to the Bulge spectrum, the 1.8 MeV line is now 
very prominent in the data (bottom- right panel). The red 
line shows the best-fitting Gaussian at 1.809 MeV with fiux 
Fi.s = (4.1 ± 0.5) 10"" phot s"^ cm"2. The blue line in the 
bottom-left figure now shows a 511 keV line with flux F^n — 
0.409 X _Fi,8 . The data do show the presence of a 511 keV line 
in the Disk component, although our experiments with the 
alternative background models and various spatial patterns 
demonstrate that the parameters of the line near 511 keV 
are not determined reliably. 
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Figure 7. Expected contribution of positrons produced by ^^AI decay to the annihilation line flux in the vicinity of the Galactic Centre. 
The light-bucket technique was used (see t|3.1ll . The slice has a width of 16° perpendicular to the Galactic plane and is centred at 6 = 0° . 
The step along the Galactic plane is 4° . The measured flux in the 508-514 keV band is shown with blue crosses; the red crosses show 
the flux in the 1804-1813 keV band scaled by a factor of 0.409, which is the expected 511 keV line flux arising from ^^Al decay under the 
assumption that 100% of positrons annihilate through the formation of positronium. It will be observed later ([[SJ that this assumption is 
supported by spectral data. Clearly the 1804-1813 kcV flux arising from the Galactic Centre region is far too low to explain the observed 
511 keV line by the local -^^Al decay. 



Thus, Fig. [S]suggests that with the two-component spa- 
tial Bulge-f Disk model: 

(i) The strong 511 keV line is confined to the Bulge 
component, with flux 0.84 - 0.93 10~^ phot cm"^ 
and the corresponding 511 keV line production rate of 
~ 0.5 10'''^ phot s~^ (equivalent to an annihilation rate of 
~ 10*'^ positrons per second if the positronium fraction is 
100%). 

(ii) The 1.8 MeV line is on the contrary confined to 
the Disk component. The flux (model flux integrated 
over a 80° x 80° square region around the GC) is 
4.1 10~^ phot s^^ cm~^, corresponding to a n ^^Al mass 
of ~ 2.8 Mq within the same region (see also IWang et al.l 
I2009D . 

(iii) An upper limit on the ratio of the 1.8 MeV and 
511 keV line fluxes in the Bulge component is -Fls/Jsh < 
0.027 {2a), which is two orders of magnitude lower than the 
expectation for the positrons being produced by '^^Al decay 
and fps ~ 1 (even if /ps — 0, the ratio is still ^^20 times 
lower than the expectation). 

(iv) The ratio of the 1.8 MeV and 511 keV line fluxes 
in the Disk component Fsn/Fi.s — 0.56 ± 0.2 is consistent 
with ^''Al decay producing the positrons, although the large 



uncertainty in the 511 keV flux in the Disk component pre- 
cludes a firmer conclusion. 



5 SPECTRA 

The emission seen by SPI from the central radian of the 
Galaxy contains contributions of many spatially and physi- 
cally different components, e.g. bright compact sources, un- 
resolved weak sources, emission coming from the disk, etc. 
The clean "decoupling" of 511 keV and 1.8 MeV emission 
into the Bulge and Disk components made in §3] strongly 
suggests that analyzing the spectrum associated with the 
Bulge component is advantageous for understanding the na- 
ture of 511 keV emission. Indeed, the Disk component might 
absorb some physically unrelated emission (like the 1.8 MeV 
line or the continuum at hundreds keV, which may be at- 
tributed to the interaction of cosmic rays with the ISM), 
and the remaining Bulge component spectrum should be 
"cleaner". The main disadvantage of this approach is that 
the significance of all features goes down once two spatial 
components with free normalizations are allowed (see e.g. 
Table[T]- the error on the 508-514 keV flux from the Bulge is 
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Figure 9. Broad 30-1000 kcV spectrum of the Bulge component of the two-component spatial model described in ^3.51 The black line 
shows the spectral model consisting of a Gaussian line at 511 keV and the associated ortho-positronium continuum, with the parameters 
given in Table O 



twice as large in the Bulge+Disk model as in the pure Bulge 
model) . In spite of the increased uncertainties, we choose to 
apply the subsequent spectral analysis to the Bulge compo- 
nent of the Bulge-|-Disk model. The broad-band spectrum of 
this component is shown in Fig. |9l In the Bulge component 
spectrum, the flux at energies above 511 keV is consistent 
with zero, suggesting that a broad band continuum does not 
contribute much around 500 keV. At energies lower than 
~300 keV, a broad band continuum is clearly visible. We 
know many compact sources that contribute to this con- 
tinuum, e.g. 1E1740.7-2942 or GRS1758-258. They have 
hard spectra around 100 keV, but decline strongly at higher 
energies, making a small contribution at ~500 keV. 



5.1 In-flight annihilation 

As a rule, the positrons produced by any plausible physical 
process are born "hot", i.e. with kinetic energy of order of 
or higher than the rest energy. For positrons with energy 
Emu ^ mef? , the rate of Coulomb and ionization energy 
losses dominates over the annihilation rate and the positrons 
first slow down before they have a chance to annihilate. 

In some scenarios (e.g. dark matter annihilation 
or production of positrons in the vicinity of relativis- 
tic compact objects), the initial energy of the positrons 
can be much higher tha n rrieC? . As discus s ed b y 



Aharonian fc AtovanI (Il981allbh : iBeacom fc Yuksell (|2006h : 
Sizun. Casse. fc Schannd (|2006l ), if -Ekin ^ m^f? , the energy 
losses and annihilation rates differ by only a factor of ~10 
and a substantial fraction of high energy positrons annihilate 



during the slow-down process. These high energy positrons 
will form a broad feature at a mean energy ~ meC^ -I- i5kin/2, 
while the slowed down positrons will annihilate at much 
lower energies and will power a narrow 511 keV line. Thus, 
a high energy (above 511 keV) component is expected to 
be present in the spectrum, with the relative intensity with 
respect to the narrow 511 keV line depending on the initial 
positron energy -Ekin and the ionization state of the inter- 
stellar medium. The latter controls the contribution of the 
ionization and Coulomb losses to the total energy loss rate. 

Shown in Fig.[TU]is the observed spectrum of the Bulge 
component and the expected in-flight annihilation spec- 
trurrjj for a neutral (solid) and ionized (dashed) medium. 
The difference in normalization is due to the larger energy 
losses in the ionized medium, which increases the fraction of 
the slowed down positrons at the expense of in-flight anni- 
hilation. The fraction of slowed down positrons annihilating 
via positronium formation was set to /ps = 0.97. The spec- 
tra shown correspond to initial positron energies of 1, 3, 
5, 10, 50 and 100 MeV. No significant flux above 511 keV 
is observed by SPI in the Bulge component (for the two- 
component model described in ii3.5p . and in Fig. [TD]we show 
the corresponding 2cr upper limits. Evidently, the SPI data 
above 511 keV do not place tight cons traints on the ini- 
tial en e rgy of positrons. As discussed by iBeacom fc Yuksell 
(|2006D : [Sizun, Casse. fc Schannc. ((200(^ ). tighter constraints 



The expected spectrum in addition to the in-flight an- 
nihilation includes the internal bremsstrahlung com ponent 
llBeacom. Bell, fc Berto"n3l2005l ; IBeacom fc Yiikse]||200^ . 
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Figure 10. Expected in-flight annihilation spectrum for a neu- 
tral (solid) and ionized (dashed) medium, normalized to the 
511 keV line flux observed in the Bulge component of the 
Bulge-|-Disk model. The initial energy of positrons is 1, 3, 5, 
10, 50 and 100 MeV. The upper limits are Icr values. Both 
the in-flight annihilation and internal brcmsstrahlung compo- 
ne nts are included in the shown spectra, calculated f ollow- 
ing I Aharonian fc AtovanI lll981bl ): [B"eacom fc Yiiksell ll2006h . The 
blue solid line shows the pure in-flight annihilation spectrum 
(without internal bremsstrahlung) for 100 MeV positrons in a 
neutral ISM. 



come from COMPTEL data combined with the SPI mea- 
surements of the 511 keV line flux and /ps, restricting i?kin 
to less than 3-7.5 M eV, depending on the ionizat ion state of 
the medium (see also I Aharonian fc Atovanll981al . for similar 
calculations based on earlier measurements of the gamma- 
ray flux from the GC region). 

The limits on the Bulge flux above 511 keV do not 
strongly constrain the amount of cosmic rays in the Bulge, 
which may produce positrons via tt^ production. Indeed, 
the same cosmic rays would also produce a comparable 
amount of tt", which would be visible as gamma- ray emis- 
sion with a peak around 100 MeV. As discussed by e.g. 
lAharonian fc AtovanI (|2000l ) , the total gamma-ray flux from 
the Inner Galaxy associated with 7r° decay does not ex- 
ceed ~ 10~^ phot cm~^ s~^ sr~^, corresponding to a flux 
< 10~^ phot cm^'^ s^^ withing the Bulge area (assuming 
a solid angle of the Bulge of ~ 0.1 sr. This is about two 
orders of magnitude smaller than the observed Bulge flux 
10"'' phot s"'^ in the annihilation line. The curves shown 
in Fig. [To] are scaled by the observed rate of positron pro- 
duction and they approximately match the SPI upper limits 
(for the initial energy of positrons ~100 MeV). Therefore, 
the SPI upper limits on the in-flight annihilation of positrons 



Table 3. Bcst-fltting parameters of the Bulge component spec- 
trum (for the Bulgc-(-Disk model) in the 400-600 kcV band. The 
quoted errors are la errors for a single parameter of interest. For 
the best-fltting model, = 401.8 for 395 degrees of freedom. 



Parameter 


Value 


Energy, keV 


510.985 ± 0.049 


Width (FWHM), keV 


2.40 ±0.16 


Flux in 511 keV line 


(0.84 ±0.03) 10-3 


Positronium fraction, /ps 


1.00 ± 0.02 



constrain the -k^ production rate to be less than ~ 10"'^ . 
Since the 7r+ and production rates by the same cosmic 
rays are comparable, it is obvious that observations of ~ 100 
MeV gamma-rays provide much tighter constraints on the 
amount of cosmic rays in the Bulge. 

5.2 Annihilation in a static ISM 

We now proceed with spectral fitting of the 511 keV line and 
ortho-positronium continuum below 511 keV. The simplest 
possible model is a combination of a Gaussian at 511 keV to 
describ e two-photon annihil ation and the three-photon spec- 
trum of lOre fc PowelJ (| 19491 ). The line normalization, energy 
and width and the normalization of the ortho-positronium 
continuum are free parameters of the model. The best-fitting 
values of these parameters are given in Table [S] The val- 
ues agree well with th e previously repor ted SPI results (e.g. 
ICh urazov et all I2OO5I : [j ean et aU I2OO6I ). Minor variations 
of the parameters among different publications can be at- 
tributed to the different background models and differences 
in the spectra extraction procedures. 

Adding a power law component does not improve the 
or affect the parameters of the model when the data in 
the 400-600 keV band are considered. As mentioned already 
at the beginning of Section [5] this is true when the Bulge 
component of the Bulge+Disk model is considered. It seems 
that the contribution of known compact sources is relatively 
small around 500 keV, while the hard continuum associated 
with cosmic-rays interactions with the gas largely goes into 

the Disk component. 

Bussard. Ramatv. fc DrachmanI 



As discuss ed by , . 

(Il979l) (see also iGuesso uni. Ra matv. fc Lingenfelteij Il99ll : 
iGuessoum. Jean, fc Gillard .2005 ). the effective line width 
and the positronium fraction principally depend on the 
temperature and the ionization state of the medium. 
This dependen ce is shown in the left panel in Fig. [TT] 
(adopted from IChurazov et al.l |2005| ). where the relation 
between the effective width of the line and the fraction of 
positrons annihilating via positronium formation is shown 
together with the observed parameters taken from Table 
[21 Compar ed to the previously reported version of the 
same plot (jChurazov et al.ll2005l ). the uncertainties on the 
positronium fraction have decreased substantially, but the 
principal result remains the same: in terms of a single-phase 
ISM model, a combination of the effective widtljj of the 



The 



effective 



vidth 



eFWHM 
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line of ~ 2.4 keV and the positronium fraction fps ~ 1 can 
be explained either with a gas of temperature in the range 
from 7500 to 4 10* K or with a much cooler gas at ~100 K. 
In both cases, the ionization degree should be at the level 
of a few %. While formally the higher temperature solution 
seems to be preferable because of the large positronium 
fraction, our experiments with background (and spatial) 
models have shown that we can drive the positronium 
fraction down to 96-97%, making the low-temperature 
solution more consistent with the data. 

Of course, the shape of the annihilation spectrum in 
an ISM of given temperature and ionization state cannot 
be fully described by a combination of a Gaussian and 
the ortho-positronium continuum. Indeed, in-flight annihi- 
lation typically contributes a broader component to the 511 
keV line than the two-photon annihilation of thermalized 
positrons. We therefore ran a grid of models for a set of 
ISM temperatures and ionization states, calculated the ex- 
pe cted annihilation spec tra using the same procedure as 
in lChurazov et aTl (|2005l ) and compared them with the ob- 
served Bulge spectrum. The initial energy of positrons was 
set to E^kin =500 keV. 

The results are insensitive to the particular value of ini- 
tial energy unless -Bkin ^ meC^ or instead -Ekin is much lower 
than a fraction of keV. Contours of values with an incre- 
ment Ax^ = 3 are shown in the right panel of Fig. [TT] . The 
minimum of x^ achieved at T = 10^ K and the ionization 
degree rj = ne/{ne+nH) = 0.025 is marked with a black dot. 
As expected, there are two pronounced minima: at ~ 100 K 
and at ~ 10* K, both at r; '-^ a few per cent. The correspond- 
ing spectra are shown in Fig. 1121 they are further discussed 
in i|5.3l below. In this analysis, the "high-temperature" is- 
land has a x^ larger than the "low-temperature" solution by 
Ax^ ~ 6. Given the simplicity of the model, this difference 
in not sufficient to make a robust conclusion in favour of the 
low-temperature solution. We view both solutions as accept- 
able. The intermediate values can be reliably excluded. This 
non-monotonic behaviour over temperature arises because 
at temperatures in the range 1000-7000 K, i) the typical en- 
ergy of a thermalized positron is below the 6.8 eV threshold 
of positronium formation via charge exchange with a hydro- 
gen atom and ii) radiative recombination with free electrons 
is suppressed relative to the direct annihilation with bound 
electron^. As a result, the positronium fraction is small at 
intermediate temperatures. 

Overall, the results of this analysis and the implications 
for the properties of the ISM broadly agree with previou s 
analyses of SPI data JChurazov et al.ll2005l:|Jean et al ]|2006l\ 
and e arlier analyses by Bussa rd. Ramatv. fc DrachmanI 



j 19791 ) :[g uessoum. Ramatv. fc Lingenfeltcn (jl99l|) 



5.3 Annihilation in a cooling ISM 

The model spectra used in the previous section were cal- 
culated assuming that the temperature and the degree of 



jGuessoum. Ramatv. fc LingenfelteJ Il99lh of the line is de- 
fined as an energy interval containing 76% of the line photons. 
For an ionization degree of a few %. 



ionization of the ISM do not change on the time scales of 
positron slow-down and annihilation. This is not an obvi- 
ous assumption since the time scale for the annihilation of 
thermalized positrons can be longer than the ISM radiative 
cooling time scale, as shown in Fig. [TS] We define an isobaric 
cooling time scale (black line in Fig. I13p as the ratio of the 
gas enthalpy to the radiative losses: 

7 nkT 
" 7 - 1 n2A(r) ' 

where 7 = 5/3 is the gas adiabatic index, n the total gas 
density, k the Boltzmann constant, T the gas temperature, 
and A(T) is t he appropriately normalized cooling function, 
adopted from [Sutherland fc Dopital (| 19931 ) . The slow-down 
time estimated for 500 keV positrons is shown with the blue 
line; for this estimate we considered only Coulomb and hy- 
drogen ionization losses (ignoring excitation losses and the 
contribution of helium) . The temperature dependence of the 
slow-down time scale in Fig. [T3] arises due to the varia- 
tions of the degree of ionization with temperature. For this 
plot we adopted the coUision al ionization equilibrium from 
lArnaud fc Rothenflud (|l985l ) . Finally, the red line shows the 
time scale for annihilation of thermalized positrons, which 
includes direct annihilation with free and bound electrons, 
charge exchange with hydrogen and radiative recombina- 
tion. 

From Fig. [13] it is clear that for temperatures in the 
range from 3 10* to ~ 3 10^ K, the gas is able to cool faster 
than thermalized positrons can annihilate. For temperatures 
outside this range, the annihilation time is shorter than the 
cooling time and the medium can to a first approximation 
be regarded as static. The most interesting conclusion that 
can be drawn from the comparison of time scales in Fig. [13] 
is that if ~ 500 keV positrons are born in a 10® K ISM, 
then the gas temperature will dro|Jf| to a few 10* K before 
the positrons have a chance to annihilate. Once the ISM has 
cooled down, the positrons will annihilate and the observer 
will measure spectral characteristics typical of annihilation 
in warm gas, with no indication that the positrons were born 
in a hot medium. 

We now proceed with calculating spectra expected for 
annihilation of positrons in a cooling ISM. This implies solv- 
ing a time dependent problem: for given initial gas temper- 
ature To and kinetic energy of positrons Eo find the annihi- 
lation spectrum emitted while the gas cools down to a given 
final temperature Ti. We solve this problem using a mod- 
ified version of our Monte Carlo code, which retains from 
the original version the algorithms for description of slow- 
ing down, thermalization and annihilation of positrons (in a 
pure hydrogen plasma). 

The outcome of this calculation crucially depends on 
the assumptions made about the ionization state of the ISM 
during its cooling from > 10® to < 10* K. In real astro- 
physical situations, a radiatively cooling plasma is expected 
to be substantially overionized relative to coUisional ioniza- 
tion equilibrium (CIE), since the characteristic cooling time 
is shorter than the time scale of hydrogen recombination. 



* In the absence of an external source of heat. 
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Figure 11. Left: The effective FWHM of t he 511 keV line versus the fraction of annihilations through positronium formation (the 
curves are adopted from lChurazov et al.|[2005l v The gray area is the width and the positronium fraction measured by SPI. There are two 
groups of theoretical curves: cold - T ^ 5000 K (dotted lines) and warm/hot — T ^ 7000 K (solid lines). The temperature is fixed for each 
curve (the labels next to the curves), but the ionization fraction varies so that plasma changes from neutral to completely ionized along 
the curve. For cool temperature curves and for the 8000 K curve, the points corresponding to the ionization degree of 0.01 and 0.1 are 
marked with the open and solid squares respectively. Each high temperature curve has two regimes, shown by thin and thick solid lines, 
respectively. Thin (thick) lines correspond to ionization fractions smaller (larger) than expected for collision dominated plasma at this 
temperature. The over-ionized state (the regime shown by thick lines) is more natural for the ISM than the under- ionized one (thin lines). 
Finally, the dashed line shows the relation between the line width and positronium fraction for a completely ionized plasma as a function 
of temperature. Right: Countours of for a grid of models over the static ISM temperature T and ionization degree 77 = ne/{ne -l-n^f). 
The initial energy of positrons is set to 500 keV. The is calculated for the 400-600 keV range for the spectrum of the Bulge component 
in the Bulge-|-Disk model. The minimum of = 402.4 achieved at T = 10^ K and r) = 0.025 is marked with a black dot. Contours are 
plotted with an increment Ax^ = 3. There are two pronounced minima: at T ~ 100 and ~ 10'* K, both with the ionization degree rj ~ 
few percent. The "high-temperature" island has a minimum x^ larger than the "low-temperature" solution by Ax^ ~ 6. This difference 
while formally statistically significant is reduced when another background model is used. Long- and short-dashed lines show the relation 
between the temperature and the degree of ionization for coUisional ionization equilibrium and radiatively cooling plasma (see i|5.3l l. 
respectively. 



This will influence both the ISM cooling rate and the fate of 
positrons. We therefore adopt the temperature dependencies 
for the non-equilibrium ionizati on state and cooling rate of 
an isq chori cally cooling gas f rom Sch mutzler fc Tscharnuten 
l|l993h and lGnat fc Sternberg ( 2007. ), which cover the tem- 
perature range from 10'^ K to 10^ K (we assume that hy- 
drogen is fully ionized at T > 5 10^ K). The temperature 
dependencies of the hydrogen ionization fraction for the CIE 
and non-CIE cases are compared in the right panel of Fig. 1111 

We carried out simulations covering a broad range of 
parameter values: Eq =10 keV-10 MeV, To = 10*-10'' K. A 
run was terminated when the ISM cooled down to 10^ K; it 
turns out that the absolute majority of positrons annihilate 
before the gas reaches this temperature, except for extreme 
values of the initial ISM temperature and positron energy. 
Although we assumed isochoric cooling, this assumption is 
not crucial, since the case of isobaric cooling is character- 



ized by very similar ISM ionizat ion state and cooling rate 
(|Schmutzler fc Tscharnuteilll993h . 

Similarly to the case of annihilation in a static medium 
considered in the previous section, we fit the spectra pre- 
dicted by our cooling ISM models to the spectrum measured 
by SPI. In Fig. 1141 we show contours of values on the [Eq, 
To) plane. The best-fitting parameters are Eq ~ 250 keV 
and To ~ 6 10^ K (see the corresponding model spectrum 
in Fig. [15] and its fit to the SPI data in Fig. [TSjl , with a 
large associated uncertainty region: the initial gas temper- 
ature and positron energy are allowed to take values from 
2 10^ to ~ 10® K and from 10 keV (the lower boundary 
of our simulations) to ~ 2 MeV, respectively. Therefore, the 
predicted gamma-ray spectrum is fairly insensitive to the 
initial parameters of the ISM and positrons. The reason for 
this is that in a broad range of Eo and To values the main 
contribution to the total spectrum is provided by annihila- 
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Figure 12. Left: Expected annihilation spectrum for a 500 keV positron in a 100 K plasma with the hydrogen degree of ionization of 0.025. 
Right: The same as in the left figure, but for 10* K plasma and 77 = 0.016. The lines show: solid black line — total annihilation spectrum; 
solid red line - total line spectrum; solid blue line — othro-positronium continuum; dotted black line - in-flight charge exchange; dotted 
red line — thermalized positron annihilation with bound electrons, dotted blue line - radiative recombination of thormalizod positrons, 
dotted green line - direct annihilation; magenta - charge exchange for thermalized positrons. 



Table 4. Values of for various spectral models for the Bulge 
component spectrum in the 400-600 keV band (400 channels, 0.5 
keV wide). 



Model 


Parameters 






Static ISM 


T = 


10^ 


K;v = 


0.025 


402.4 


Static ISM 


T = 


10* 


K; r,= 


0.016 


408.6 


Cooling ISM 


To = 


= 6.3 


105 K; 


Eo = 250 keV 


400.8 



tions of thermalized positrons during the gas cooling from 
~ 10^ K to ~ 10* K (when the positron annihilation time is 
a few times the ISM cooHng time, see Fig. I13p . 

Judging by the small differences between the minimum 
values, summarized in Table 3] (and recalling again that 
there are non-negligible systematic uncertainties associated 
with the data analysis), our best-fitting static ISM and cool- 
ing ISM models describe the SPI spectrum nearly equally 
well. In other words, we cannot strongly favour any of these 
models based on the available spectral data. In reality, this 
is mostly owing to the finite energy resolution of SPI. This 
becomes clear if we compare the spectra predicted by the 
different models, shown in Figs. I12land ll5l 



6 DISCUSSION 

A comprehensive answer to the fundamental question of the 
origin of the bulk of annihilating positrons should explain a 
number of key observational facts, including the total pro- 
duction rate of positrons (~ 10*^ positrons per second within 
the Galactic bulge), the dominance of the bulge compared 
to the disk, the morphology of the 511 keV line emission 
(a 8° X 6° "Bulge" plus a more extended, along the Galac- 
tic plane, "Disk" component) and the spectral properties of 
the Bulge component, which (within the accuracy allowed 
by the SPI energy resolution) resemble the annihilation of 
positrons in a warm (~ 10'' K) or cold (~ 10^ K) medium 
with the ionization degree of a few per cent. 

On these grounds (flux and morphology) , such straight- 
forward explanation as ^^Al or the interaction of cosmic rays 
with the ISM are disfavored. The morphology problem can 
be partly alleviated by assuming transport of the positrons 
away from the place where they are born (e.g. IPrantzod 
l2006h . Recently, an at tempt to provide a self-consistent 
model was done by [Higdon. Lingenfelter. fc Rothschildl 
(|2009l ): iLingenfelter. Higdon. fc Rothschildl (|2009l ). Their 
model assumes that supernovae are the source of positrons 
(via 13+ decay of ^^Ni, ''*Ti and ^^Al) in the Galaxy and the 
propagation of positrons in the ISM is different in the bulge 
and in the disk. The difference in propagation explains the 
weakness of the disk compared to the bulge - a large fraction 
of positrons born in the disk escape into the halo, while only 
a small fraction of positrons escape the bulge. In the bulge. 



© 2001 RAS, MNRAS OOO.[Tlf20l 



Positron annihilation spectrum from the Galactic Centre 17 





1 1 ' 


1 ' ' ' ' 








: o 






lati 












1 s 












i < 








Cooling 

1 1 I 


Annihilation : 



10* 10^ 10« W 10" 



T. [K] 



Figure 13. Comparison of different time scales relevant for an- 
nihilation of positrons in a cooling ISM as a function of its tem- 
perature: cooling of gas (black) , slowing down of positrons (blue) , 
annihilation of thermalized positrons (red). Calculations are done 
for the proton density 1 cm~'^. The initial energy of positrons is 
assumed to be 500 keV. The dotted vertical lines separate the 
areas where cooling or annihilation dominates. For temperatures 
in the range from ~ 3 10* to ~ 3 10^ K, the gas can cool faster 
than thermalized positrons annihilate. For instance, if 500 kcV 
positrons are born in a ~ 10® K ISM, the gas will cool down 
to a few 10* K before the positrons have a chance to annihilate. 
Once the temperature has dropped, the thermalized positrons will 
annihilate and the observer will see a gamma-ray spectrum indi- 
cating annihilation in a ~ 10* K gas, regardless of the initial ISM 
temperature. 



the positrons diflFuse through the ISM and annihilate when 
entering the HII and HI envelopes of molecular clouds - this 
explains the spectral properties of the annihilation emission. 
Finally, the morphology of the bulge is explained by anni- 
hilation in a Tilted Disk within 1.5 kpc of the GC. The key 
assumptions in the Higdon et al. model are that i) mag- 
netic flux tubes are almost radial in the Bulge region and ii) 
positrons propagation along the flux tubes can be treated in 
the diffusion approximation with the frequency of scattering 
(or the mean free path) estimated from scaling the energy 
density of magnetic fields down to the Larmor radius. Their 
estimate of the mean free path in the Bulge region is large 
enough to allow positrons to reach the Tilted Disk (see tj3.4p 
before they slow down and annihilate. Thus the morphology 
of the Tilted disk gives the source of annihilation radiation 
in the GC its characteristic shape. 

Here we consider the alternative possibility that the 
thermal evolution of the ISM plays the key role in the GG 
annihilation emission. The key assumption of the model is 




10" 10^ 10^ 

Eo, eV 

Figure 14. Countours of for a grid of models for a cooling ISM 
over the range of initial ISM temperatures Tq and initial positron 
energy Eq. The minimum of = 400.8 achieved at Tq = 6.3 10^ 
K and Eo = 250 keV is marked with a black dot. Countours are 
plotted with an increment Ax^ = 3. 



that the medium in which the bulk of positrons are produced 
and annihilate has a cooling time shorter than the slow down 
and annihilation time scales. From Fig.[l3]it follows that the 
temperature ought to be lower than 3 lO^K. 

Gas with temperature in the range 10^- few 10® K is 
observed in the Galaxy, as indicated by e.g. CVI, OVII and 
OVIH lines fe.g. lMcGammon et all2002l ). An ISM with tem- 
perat ure ~ 3 10® K is also obse rved in the central region of 
M31 (|Bogdan fc Gilfanovllioosi ) . The origin of this gas is not 
clear. For a long time observations of 6.7 and 6.9 keV lines 
of heavily ionized iron and an unresolved X-ray glow (at en- 
ergies above a few keV) of the Galactic plane were regarded 
as evidence for a much hott er, ~ 10^ K, plasma filling most 
of the bulge's volume (e.g. iKovama et al.lll98d ). Positrons 
injected into such a medium will slow down and annihilate 
before the plasma is able to cool (Fig. I13|l , and a very broad 
annihilation line without an ortho-positronium continuum 
is expected to be observed, in stark contrast with observa- 
tions. However, recent observations (iRevnivtsev et al.l l2006l : 
ISazonov et al]l2006l : IRevnivtsev et al.ll2009^ have convinc- 
ingly demonstrated that the apparently diffuse Galactic X- 
ray emission is in fact a superposition of millions of faint 
compact sources - accreting white dwarfs and coronally ac- 
tive stars. This removes the problem of confining ~ 10* K 
gas by the gravitational well of the Galaxy. Cooler, 10^- few 
10® K, gas can on the contrary be trapped in the Galaxy 
potential. For instance, the density of 2 10® K gas in hydro- 
static equilibrium in the Milky Way potential would drop 
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Figure 15. Left: Expected annihilation spectrum for a radiatively cooling plasma with initial temperature To = 6.3 10'' K and initial 
positron kinetic energy Eq = 250 keV. Different lines show different spectral components, as in Fig. 1121 In contrast to the annihilation in 
a static medium, in-flight component makes negligible contribution for the cooling ISM solution, in which case the annihilation spectrum 
is essentially a superposition of spectra generated by thermalized positrons in gases with temperatures ranging from ~ 10^ to ~ 2 10^ K. 
Right: Spectrum of the annihilation radiation for a cooling ISM with initial temperature To = 6.3 10^ K and initial positron energy 
Eq = 250 keV in comparison with the observed Bulge component spectrum in the Bulge+Disk model. 



by a factor of ~ 3 10^ at a distance of ~ 1.5 kpc above 
the plane. This means that such gas can accumulate in the 
Galactic bulge. 

However, a self- consistent model explaining the few 
10** K ISM in the bulge of the Milky Way or in M31 is still 
miss ing. If one takes the mass loss rate by evolving stars 
(e.g. ICiotti et al.lll99l|) and the total ener gy and iron input 
rate from SNIa (e.g. Mannucci et al.|[2008l ). then a "natural" 
temperature of the gas is of order few keV (mean energy 
per injected particle) and the expected abundance of iron 
is about 5 times the Solar value. These parameters do not 
fit the p roperties of the ob serv ed few 10^ K compon ent of 
the ISM. lTang et all (|2009l ') and lTang fc Wand (|20ld ) argue 
that the account for non-uniformity and intermittence in en- 
ergy injection by SNIa leads to a broad gas distribution over 
temperature with a lower than the mean temperature com- 
ponent making largest contribution to the observed X-ray 
emission. It is not yet clear if this provides a full explana- 
tion of the origin of 10® K gas. We therefore make an ad hoc 
assumption that a few million degree medium is widespread 
in the Milky Way bulge. 

The positron annihilation in a cooling ISM ( t|5.3p can 
resolve some of the issues mentioned at the beginning of 
this section. Indeed, if the ISM is able to cool down to 
10^-10* K before the positrons annihilate, then the observer 
will always see a spectrum characteristic of annihilation in 
a warm plasma. This resolves the question of why positrons 



are apparently annihilating in a warm/cold medium even 
though such a medium should have a relatively small filling 
factor. Otherwise, the positrons would have to be produced 
preferentially in the warm and/or cold ISM phases or to 
migrate between phases. In the cooling ISM model this is 
not required. If the ISM is dynamic and goes through cool- 
ing/heating cycles, then even those positrons born in the 
hot phase will annihilate in a warm (10''-10^ K) gas. Trans- 
port of positrons over large distances is not required in this 
model. Some moderate transport of positrons will not con- 
tradict the model, since it will just make shorter the lifetimes 
of positrons initially born in the hot ISM, but the positrons 
will still predominantly annihilate while migrating through 
the warm/cold ISM phases. 

Putting all arguments together, it seems plausible that 
the following model might work: 

• Positrons are injected in volume-weighted fashion into 
the ISM (e.g. by SNIa). This means that most of the 
positrons are initially in the hot ISM phase. The key re- 
quirement of the model is that the large volume fraction has 
to be occupied by a few 10*^ K gas - which may not be true 
in the Galactic bulge. 

• The plasma cools radiatively before the positrons slow 
down. An outflow, adiabatic expansion of the gas ad- 
vect/move positrons and may contribute to positron cool- 
ing. 

• Eventually, the plasma reaches relatively low temper- 
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atures and the positrons annihilate, producing the spectral 
features as observed by SPI. 

If correct, this model predicts that there should be some 
similarity in the spatial distribution of the annihilation emis- 
sion and the thermal emission of plasma with temperature in 
the 10'*-2 10® K range. Outflows from the central regions of 
normal galaxies are well-known. A good example is M31 with 
clear evidence of hot ~ 3 10*^ K gas (|Bogdan fc Gilfanovl 
I2OO8I ). possibly outflowing from the central region. In sub- 
sequent publications we plan to consider the dynamics of 
galactic outflows, advection and annihilation of positrons 
and the expected spectral signatures (in soft X-ray and other 
wavebands) of the main annihilation sites. 



7 CONCLUSIONS 

We have analysed SPl/INTEGRAL data on the 511 keV 
line from the Galactic Centre, accumulated over ~6 years of 
observations. 

Using a simple two-component Bulge-fDisk model we 
have constructed the broad-band energy spectra of both 
components. The Bulge component shows a prominent 511 
keV line and essentially no flux at 1.8 MeV, while the Disk 
component on the contrary contains a prominent 1.8 MeV 
line and a very weak annihilation line. The inferred annihi- 
lation rate in the Bulge component is ~ 10*^ positrons per 
second (for fps = 1). 

The FWHM of the 511 keV in the Bulge is 2.40 ± 0.16 
keV. The flux ratio of the line and three-photon contin- 
uum corresponds to a fraction of positrons annihilating via 
positronium formation close to unity (/ps = 1.00±0.02). Our 
experiments with using slightly different SPI background 
models have yielded fps in the range 0.96-1. This combi- 
nation of spectral parameters can be explained in terms of 
a single-phase static ISM model, namely as annihilation of 
positrons in a warm (~ 10* K) or cold 10^ K) medium 
with the hydrogen ionization fraction of a few per cent. 

We have demonstrated that the observed annihilation 
spectrum of the Bulge component can also be explained with 
a model of positrons annihilating in a radiatively cooling 
ISM, provided that the initial temperature of the ISM is in 
the range 2 10^-10^ K. In such a medium, the ISM cools 
down to ~ 10^-10* K before the positrons are able to anni- 
hilate and the observer sees an annihilation spectrum typi- 
cal of a warm ISM. This suggests that the positrons might 
be initially injected into the hot ISM phase, which has the 
largest filling factor in the Galaxy, but the annihilation oc- 
curs only when the gas has cooled down to a relatively low 
temperature, effectively erasing the memory on the initial 
conditions (the initial temperature and ionization state of 
the ISM and the initial energy of positrons). If there is an 
outfiow of hot gas, then the positrons will be advected with 
the flow and annihilate in the part of the flow where the 
temperature has dropped to 10^-10* K. The mass of the 
cooling gas need not to be large - it is enough to assume 
that the initially hot ISM phase occupies the largest fraction 
of the volume where the positrons are produced, provided 
that positrons do not migrate freely between phases. The 



key requirement of the model is that the large fraction of 
the bulge's volume is occupied by a few 10^ K gas. In this 
paper this is adopted as an ad hoc assumption. The origin 
of this gas in not entirely clear, although observations sug- 
gest that this component is present in the central part of the 
Galaxy. 
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